SWIMMING AS A LIMIT CYCLE 

HENRY O. JACOBS 



Abstract. Steady swimming can be characterized as both periodic and stable. 
These characteristics are the very definition of limit cycles, and so we ask "Can 
we view swimming as a limit cycle?" In this paper we will find that the answer 
is "yes". We will define a class of dissipative systems which correspond to 
the passive dynamics of a body immersed in a Navier-Stokes fluid (i.e. the 
dynamics of a dead fish). Upon performing reduction by symmetry we will 
find a hyperbolically stable fixed point which corresponds to the stability of a 
dead fish in stagnant water. Given a periodic force on the shape of the body 
we will invoke the persistence theorem to assert the existence of a loop which 
approximately satisfies the exact equations of motion. If we lift this loop with 
a phase reconstruction formula we will find that the lifted loops are not loops, 
but stable trajectories which represent regular periodic motion reminiscent of 
swimming. 
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1. Introduction 

The motion of steady swimming appears periodic [Sha98]. Thus the question: 

Can we reasonably interpret steady swimming as a limit cycle? 

In this paper we will assert that the answer to this question is "yes". But why 
would we want to interpret swimming as a limit cycle? Well firstly, a limit-cycle in- 
terpretation builds upon an existing body of knowledge derived from experimental 
and computational observations |AS05j . [L BLT03a] . Moreover, this interpreta- 
tion could be of interest to control engineers because robust and regular behavior 
reduces the complexity of controllers. 
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Figure 1. The periodic motion of a jellyfish. Photo taken from 
|KCDClT] courtesy of Kakani Katija Young. 

Main Contributions. We will understand the system consisting of a body im- 
mersed in a fluid as a dissipative system evolving on a space A. One observes the 
system is invariant with respect to the group of rigid rotations and translations, 
SE(g?). This observation suggests that one can describe the system evolving on the 
quotient manifold gw^y- The main contributions of this paper are: 

• We prove the existence of a hyperbolic stable point in the quotient space 



SE(d) ' 

• We then add a periodic potential energy which models the periodic con- 
traction of muscles. If A were finite dimensional, then adding this periodic 
perturbation would transform a hyperbolic stable point into a stable limit 
cycle. However, A is not finite dimensional. Nonetheless, we will prove 
the existence of loops in g^y which approximately satisfy the equations of 
motion to arbitrarily good accuracy. 

• We prove that the dynamics in A which corresponds to loops in 

are given by regular periodic motions, where each period is related to the 
previous by a rigid rotation and translation. 

This all suggests that orderly behavior such as steady swimming in periodically 
perturbed systems could be the rule and not the exception. 

In greater detail, the paper proceeds as follows. We will describe the motion 
of a flexible body immersed in a Navier-Stokes fluid as a dissipative system on a 
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tangent bundle TQ. This passive, or unactuated system, will be referred to as the 
"dead fish" system to specialize the terminology to the situation at hand. We will 
then find this systems exhibits two symmetries. The first symmetry corresponds to 
the invariance with respect to the particle-labeling of the fluid and is represented 
by a Lie group G. Reduction by G will send us from dynamics on TQ to dynamics 
on (TQ)/G. The second symmetry is invariance with respect to global rotations 
and translations of the system, and will send us from dynamics on (TQ)/G to 
SE^fT • O n this final space we will find a hyperbolic stable point. We will then 
add a periodic force which acts on the shape of the body to model periodic muscle 
contractions of a fish. Then we will attempt to use the persistence theorem to 
obtain a stable limit cycle in this periodically perturbed system. I say "attempt" 
because we will fail. To get past this failure we will be aim for a more modest 
goal and successfully use persistence theorem to assert the existence of a loop 
which approximately satisfies the exact equations of motion (to an arbitrarily good 
accuracy). Upon lifting a given loop in to (TQ)/G we will find regular 

motion whereupon each period of the perturbation is a constant rigid rotation 
and translation of the previous period. Such behavior is reminiscent of steady 
swimming. Schematically, this idea is represented via the commutative diagram 



i i r> i periodic pertubation 

dead fash > swimming 



(1) *SE(3)\ 



n SE(3)\ 



dead fish periodic perturbation perturbed dead fish 
SE(3) y SE(3) 

The top left corner of diagram ([I]) represents a dead fish in water as a dissipative 
mechanical system. If we reduce the system by SE(3) we will get a system on a 
quotient space symbolized by the bottom left corner of diagram (|lp. In this bottom 
left corner we find a hyperbolic stable point which asserts "a motionless body in 
stagnant water is a stable state for dead fish". The horizontal arrows of diagram 
(jl| correspond to the addition of a time periodic pertubation on the shape of the 
body. This perturbation produces a stable limit cycle in the bottom right corner 
of diagram via the persistence theorem^ Finally, the existence of this periodic 
orbit in the bottom right corner of diagram ([I]) implies the existence of a family of 
curves which one could reasonably refer to as "swimming" in the top right corner 

ofQ. 

1.1. Background Material. There presently exists a substantial body of knowl- 
edge in the form of computational and biological experiments which are consistent 
with and support the hypothesis that swimming could be interpreted as a limit 
cycle. For example, experiments involving tethered unactuated (i.e. dead) fish 



1 There are analytic concerns however, and we deal with them in ijZ] 
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immersed in a flow behind a bluff body suggest an ability to passively harvest 
energy from the surrounding vorticity of the flow. The same studies also provide a 
relevant example of oscillatory behavior as a stable state for an unactuated system 
(e.g. that of a dead fish) |LBLT03aj . |LBLT03b| . 

For the case of live fish, periodic motor neuron actuation has been recorded di- 
rectly and periodic internal elastic forces have been approximated via linear elas- 
ticity models |Sha98] . Additionally, numerical experiments involving rigid bodies 
with oscillating forces suggest that uniform motion (i.e. flapping flight) is an 
attracting state for certain pairs of frequencies and Reynolds numbers [AS05] . Fi- 
nally, experiments involving flexible paper in vertical oscillating flows suggest that 
the vortices shed have the apparent result of stabilizing a top-heavy body despite 
steady state analysis which would suggest instability |LRW + 12j . This body of 
knowledge suggests further investigation into the role of non-stationary flows in 
steady swimming via undulatory motion. In particular, we will be investigating 
non-stationary flows which result from periodic perturbations of muscle fibers. 

In order for us to study non-stationary flows we must have a language which 
is efficient so that our simple idea is expressed with equal simplicity. In addition, 
this language must be rich enough to express fluid-solid systems. We find that the 
language of gauge theory is a particular good choice. Specifically, we will invoke 
the standard bundle structure of guage theory by setting the gauge symmetry 
equal to the particle relabeling symmetry of the fluid. This was first investigated 
in detail for the case of propulsion in Stokes fluids |SW89] and potential flow 
in ideal fluids |Kel98j and |KM00j . In these gauge theories, the shape of the 
body was controlled explicitly and the locomotion which resulted was seen as a 
consequence of the curvature of a physically meaningful principal connection. The 
framework of [Kel98j for arbitrary bodies was specialized to articulated bodies in 
[KMRMH05] and used for motion planning in |MRR06j . More recently. JVKM09] 
was able to extend [KMOOJ to flow with point vortices with the use of Lagrange- 
Poincare reduction as performed in |CMR01] . Finally, recent work has expanded 
this picture to Euler and Navier-Stokes fluids with arbitrary circulation through 
the use of Lie groupoids |JV12j in analogy with the classic work of V I Arnold on 
Lie groups |Arn66] . 

1.2. Organization and Motivation of our approach. The objective of this 
paper is to demonstrate that swimming can be interpreted as a limit cycle. Specif- 
ically, we will show that this interpretation may be viewed as a periodic pertuba- 
tion of the stable state corresponding to a dead fish in stagnant water. Therefore, 
we will begin this paper with a discussion in £j2] on limit cycles and demonstrate 
how a periodically perturbed hyperbolic stable points yields a limit cycle. Because 
the first step in this journey is the identification of a hyperbolic stable point, we 
will seek such a stable point in the realm of fluid-solid systems. This search for 
stable points begins with the Lagrangian mechanical description of fluids (jj|3.1 ) 
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and solids ( |3.2 ) and finally fluid-solid interaction (S |3.3[ ). However, despite hav- 
ing the Lagrangian formalism specified we will not find a hyperbolic stable point 
until we have made further preparations. Our system is too large and there are a 
number of redundancies which must be removed. The two redundancies which will 
be "quotiented away" are the particle labels of the fluid (£3.3) and the frame of 
reference (£3.4). Finally, our quest for a hyperbolic stable point comes to an end 
in section £ 4j where we identify one in the quotient space. With this fixed point 
we will be able to use the "persistence theorem" to assert the existence of loops 
which approximate the exact dynamics on the quotient space. Finally, in £j5] we 
will apply a phase reconstruction formula lift these loops to unquotiented spaces. 
The lifted loops will no longer be loops, but stable trajectories which represent 
regular periodic motion, generally reminiscent of swimming. 



1.3. Conventions and Notation. All objects and morphisms will be assumed to 
be sufficiently smooth. Moreover, we will not address the existence or uniqueness 
of solutions for fluid structure systems. In particular we will assume that solutions 
exist for all time. If M is a smooth manifold then we will denote the tangent 
bundle by tm '■ TM — > M. The set of vector fields on M will be denoted by X(M). 
More generally, given a fiber-bunde it : E M, we denote the set of sections of E 
by T(E). We will denote the unit circle by S 1 and the Lie group of rigid rotations 
and translations of W 1 by SE(cf). 

1.4. Acknowledgements. The italicized question in the second line of the in- 
troduction was posed to me by my wife, Erica J. Kim, while she was studying 
humming birds. Additionally, Sam Burden, Ram Vasudevan, and Humberto Gon- 
zales provided much insight into how to frame this work for engineers. I would 
also like to thank Professor Shankar Sastry for allowing me to stay in his lab for 
a year and meet these people who are outside of my normal research circle. At 
the same time, the original version of this paper was written in the context of Lie 
groupoid theory, where the guidance of Alan Weinstein was invaluable. Much of 
the groupoid theory is hidden in the present publication, and I hope to reveal the 
groupoid theory explicitly at a later date. Jaap Eldering and Joris Vankerschaver 
have given me more patience than I may deserve by reading my papers and check- 
ing my claims. The final draft of this paper was solidified with the help of Darryl 
Holm by multiple long meetings going through the paper one line at a time. This 
research has been supported by the European Research Council Advanced Grant 
267382 FCCA and NSF grant CCF-1011944. 



2. Limit Cycles 

Let M be a finite dimensional Banach manifold and let X e X(M). Moreover, 
assume we have found a hyperbolically stable fixed point xq G M. Then we 
may embed this autonomous dynamical system in the time periodic augmented 
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phase space S 1 x M with the coordinates (t, x) by using the vector field (d t , X) G 
X(S' 1 x M). If we do this then the loop 

(2) r = {(t,x ) : t G S 1 } C S 1 x M 

is a hyperbolically stable limit cycle of (d t ,X) (see figure [2] B). In fact, T is a 





Figure 2. (A) Depicts a one dimensional system with a hyperbolic 
stable point. The embedding of the system in (A) into time-periodic 
augmented phase space exhibits a limit cycle as shown in (B). Fi- 
nally, a limit cycle will persist under a time-periodic perturbation 
as a diffeomorphism of the original limit cycle. This is drawn in (C). 

special case of a normally hyperbolically invariant submanifold. Therefore under 
sufficiently small perturbations (d t ,X) i— > (d t ,X + SX) there will still exist a hy- 
perbolically stable limit cycle, T e , in the perturbed system which is close to r . 
This is a consequence of the persistence theorem |Fen72l IHPS77] . 

It is notable that the first step in going down this path is to find a hyperbolic 
stable point. We will do this for the case of a dead fish as follows. In the next 
section we will derive the equations of motion for a dead fish as an instance of 
the Lagrange-d'Alembert equations on a tangent bundle TQ. We will then find 
two symmetries of the system. One corresponding to a Lie group G and another 
corresponding to the Lie group SE(d). Upon performing reduction by symmetry 
we will obtain dynamics on the quotient space 

TQ/G 



[TQ/G] :- 



SE(d) 



we will find that the vector field on [TQ/G] has a (weak) hyperbolic stable point, 
x G [TQ/G]. We could then construct a (weak) hyperbolic stable limit cycle given 
by equation Q with M = [TQ/G]. 
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Unfortunately, Q is an infinite dimensional Frechet manifold. Therefore the 
persistence theorem does not directly apply However, we will be able to assert 
the existence of a loop which approximates the exact dynamics. Using a phase 
reconstruction formula we will be able to solve for the motion of the fish. 

3. Passive Dynamics 

In this section we construct a Lagrangian formalism for fluid structure interac- 
tion. Recall that given a Lagrangian, L : TQ — > R, the equations of motion of the 
corresponding Lagrangian system are given in coordinates by the Euler- Lagrange 
equations 

D fdL\ dL 

Dt \d^) ~ ~dq ~ ' 
Hamilton's principle then state that a solution curve, q : [0,T] — > Q, between 
% — ?(0) and qi = q(T) satisfies the Euler-Lagrange equations if and only if q(-) 
extremizes the action 

T 



%(•)]= / L(q,q)dt 
Jo 



with respect to variations of the curve q(-) with fixed end points. In otherwords, 
5S = with respect to variations 5q(t) with 5q(0) = and 5q(T) = 0. More gen- 
erally, we may consider a force field, F : TQ — > T*Q, and consider the Lagrange- 
d'Alembert equations 

D fdL\ dL 



. , F - 

Dt \dq J dq 

which are equivalent to the Lagrange-d'Alembert variational principle 

5S= [ (F(q,q),Sq)dt 
Jo 

with respect to variations Sq with fixed end points [MR99, Chapter 7]. If Q is 
equipped with a Riemanian metric, (■, -)q : TQ © TQ — > IR, then it is customary 
to consider Lagrangians of the form 

L(Q,Q) = \(qA)Q-U{q) 

where U : Q — > M. We call such a Lagrangian a kinetic minus potential Lagrangian. 
The Euler-Lagrange equations in this case will take the form 

Where is the Levi-Cevita covariant derivative and V£7 is the gradient of U 
jAMOOl Proposition 3.7.4]. Moreover, given the force F : TQ -> T*Q, the 
Lagrange-D'Alembert equations take the form 

^ t =VU(q) + UF(q,q)) 
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where Jj : T*Q — > TQ is the sharp operator induced by the metric. In the following 
subsections we will illustrate how to understand a body immersed in a fluid as 
a kinetic-minus-potential Lagrangian system with a dissipative force field in the 
above sense. 

3.1. Fluids. Consider the manifold IR 3 with the standard flat metric and coor- 
dinates x,y,z. The flat metric induces the volume form rfvol = dx A dy A dz. 
One can consider the infinite dimensional Lie group of volume preserving dif- 
feomorphisms, Diff vo i(IR 3 ), where the group multiplication is simply the com- 
position of diffeomorphisms. The configuration of a fluid flowing on IR 3 rela- 
tive to some reference configuration is described by an element <p G Diff vo i(IR 3 ). 
Given a curve <p t G Diff vo i(IR 3 ) one can differentiate it to get a tangent vector 
dp t = ^ip t G T Diff vo i(IR 3 ). One can interpret (p as a map from IR 3 to TIR 3 by the 
natural definition ip(x) = ^<p> t (x). Therefore, a tangent vector, (p G T Diff vol (IR 3 ), 
over a diffeomorphism (p G Diff vo i(IR 3 ) is nothing but a smooth map <p : IR 3 — > TIR 3 
such that Tr3 o tp = tp where r K 3 : TIR 3 — > IR 3 is the tangent bundle projection. 
This implies that (p o ip' 1 is a smooth divergence free vector field on IR 3 . We call 
(p the material representation of the velocity, while (p o ip" 1 g X vo i(IR 3 ) is the spa- 
tial representation (see figure [3]). The Lagrangian, L : T(Diff vo i(IR 3 )) — > IR, is the 
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Figure 3. A cartoon comparing the material vs spatial represen- 
tation of fluids 
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kinetic energy of the fluid, 



L{<p,<p):=± I 11^)11^ vol. 
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M 



One can derive the Euler-Lagrange equations on Diff vo i(IR 3 ) with respect to the 
Lagrangian L to obtain the equations of motion for an ideal fluid. However, this 
Lagrangian exhibits a symmetry. 

Proposition 1 QArn66]). The Lagrangian L is symmetric with respect to the right 
action Diff vol (M 3 ) on TDiff vol (M 3 ). 

Proof. Let ip G TDiff vo i(IR 3 ) over the diffeomorphism ip G Diff vo i(IR 3 ). Morever let 
if) G Diff vo i(M 3 ). Then the map ip o ip is a vector in T Diff vo i(IR 3 ) over the point 
tpoip. This defines a right Lie group action of Diff vo i(IR 3 ) on T Diff vo i(IR 3 ). We find 

L(ip o ip) = - j ||0 o ip(x) \\ 2 d vol 
2 J M 

_ 1 

~ 2 

_ 1 

~ 2 

_ 1 

~ 2 

= L{<p) 



x)||>*dvol 
||V90r)|| 2 det(</>)dvol 

y(x)\\ 2 dvo\ 

M 



which illustrates that L is symmetric with respect to the natural (right) action of 
Diff vol (M 3 ) on TDiff vol (M 3 ). □ 



Noether's theorem asserts the existence of a conserved quantity as a result of 
the symmetry of L. This conserved quantity is the circulation of the fluid and 
is known as Kelvin's circulation theorem. Additionally, this symmetry suggests 
we can write equations of motion on the quotient space TDiff vo i(M 3 )/Diff vo i(IR 3 ) 
which one can identify with the set of divergence free vector fields, 3£ vo i(lR 3 ). It 
was the discover of [Arn66] that these equations could be written as 

du 

— — h u ■ Vu = — Vp , div(w) = 0. 
at 

One should recognize as these as the inviscid fluid equations. Moreover, if we 
define the linear map, : 3£ vo i(IR 3 ) — > £*oi( 



(/^(w), w) = (i / (Au,w)dvo\ 



M 



then we derive the Lagrange-D'Alemnbert equations by lifting / to obtain a force 
field F : T(Diff vol (M 3 )) -»■ T*(Diff vol (M 3 )). If we do this, then u{t) satisfies the 
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Navier-Stokes equations 

— — h u • Vu = —Vp — fiAu , div(tt) = 0. 

Of course there is no reason that we must restrict ourselves to the case of fluids 
on IR 3 . We can do virtually the same procedure on IR 2 or any finite dimensional 
Riemannian manifold for that matter. We refer to [AK92] for details on these 
generalizations. 

3.2. Solids. Let B be a compact manifold with boundary dB and volume form 
d volg. Let Emb(£>) denote the set of embeddings of B into M. d . If Emb(i3) is non- 
empty then it is an infinite dimensional manifold which serves as the configuration 
manifold for the theory of elasticity. For example, the theory of linear elasticity 
assumes B to be a Riemannian manifold with metric, (■, -)g and uses the potential 
energy 

linear (6) = \ I trace ([/ - C b ] T ■ [I - C 6 ]) Mvol B . 

1 Jb{B) 

where C& is the right Cauchy-Green strain tensor |MH83| . This is just an example 
and we will not need our potential energy to be of this form, but this would hold 
the properties required for the main theorem of our paper. For our purposes, the 
most important property of U is (left) SE(c?) invariance and an isolated minima in 
shape space with a nondegenerate second variation. Let us first discuss what we 
mean by SE(g() invariance. Note that the standard action of SE(c?) on IR d induces 
a left action on Emb(B) by composition of maps. Specifically, given azG SE(d) 
and a b G Emb(£>) we define the embedding z ■ h G Emb(£>) by 

(z ■ b)(x) = z ■ (b(x)) ,\/xeB. 

To say that U is SE(d) invariant means U(z ■ b) = U(b) for all z 6 SE(3). This 
implies the existence of a function [U] : [Emb(S)] — > IR where [Emb(5)] := 
SE(d)\ Emb(_B) is the quotient space, also referred to as the shape space. We will 
assume that [U] has an isolated minimizing shape s m i n = argmin([£/]) G [Emb(i?)] 
and that the second variation 5 2 [U] is nondegenerate at s min . Note that s min is si- 
multaneously a point in [Emb(S)] and a submanifold of Emb(B). The Lagrangian 
is given by 

L B (b,b) := \ [ p B (x)\\b(x)\\ 2 dvo\ B (x)-U(b). 

From Lb one may compute the equations of motion for an elastic body moving 
in a vacuum. Such a system would be conservative, and perhaps unrealistic. To 
amend this we will include a dissipative force, which is given by a vector bundle 
map, Fb : T[Emb(B)] —> T*[Emb(£>)], such that the tensor 

(Ffl(-), •) : T[Emb(£)] © T[Emb(B)} -> R 

is a (weakly) positive definite form on [Emb(£>)]. Such a force has the effect of 
dampening the rate of change in the shape of the body but it will not dampen 
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motions induced by the action of SE(d). In other words, we assume that a jiggling 
body eventually comes to rest at a shape s m i n by the dissipation of energy 

3.3. Fluid-solid interaction. Let £>, Emb(£>), Lb, F& be as described in the 
previous section. Given an embedding b G Emb(i3) let rr^ denote the set 

rr^ = closure 

The set is the region which is occupied by the fluid given the embedding of 
the body b. If the body configuration is given by 6 £ Emb(S) at time t = and 
b\ G Emb(i3) at time t — 1, then the configuration of the fluid is given by a vol- 
ume preserving morphism from rr^ to rr^, i.e. an element of Diff vo i (RRi& , RRi&J. 
Given a reference configuration &o G Emb(i3) for the body we define the configu- 
ration manifold 

Q:={(b,(p) | &GEmb(£), 

if G Diff vo i (RRi 6o ,RRi fe ) 

lim (<f(x) — x) = 0}. 

x— >oo 

One should not that the manifold Q has some extra structure. In particular, the 
Lie group G := Diff vo i(RR!fe ), represents the symmetry group for the set of particle 
labels and acts on Q on the right by sending 

(b,<p) G Q \-> (b,ipoijj) e Q 

for each ip G G and (6, cp) G Q- 

Proposition 2. The projection (3 : Q — )■ Emb(i3) defined by f3(b,ip) = b makes Q 
into a principal G-bundle over Emb(B). 

Proof. The C7-orbit of (6, 99) is the set 

{6} X Diff vol (RRi bo ,RR^) C Q. 
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Therefore, by identifying each G-orbit, {b} x Diff vo i(f«y& , rs«&), with b G Emb(£>) 
it is clear that the set of G orbits is isomorphic to the set Emb(£>). Thus we find 
that Q/G = Emb(B) and the quotient map is given by /3(b, if) — b. □ 

Now we must define the Lagrangian. To do this, it is useful to note that the 
system should be invariant with respect to particle relablings of the fluid and so 
the Lagrangian should be invariant with respect to the right action of G on TQ 
given by 

(b, b, ip, ip) G TQ h-> (b, b, ip o ip, ip o ip) e TQ 

for each ip G G. As a result we can define a Lagrangian on the quotient space 
TQ/G. Incidentally, this quotient space is much closer to the space typically 
encountered in fluid-structure interaction because the fluid is almost represented 
by a vector field rather than a diffeomorphism. We say "almost" because the fluid 
velocity at the boundary is generally not tangent to the fluid domain and so this 
is technically not a vector field. To be precise, the fluid is represented by a smooth 
map from the fluid domain to vectors above the fluid domain, i.e. an element of 

r(T^ 6 R 3 ). 

Proposition 3. The quotient space TQ/G is the set 

TQ/G := {(b,b\u) | (b, b) G TEmb(B), 

u G T(Tp~ b R d ), 

(3) u(b(x)) = b(x)Vx G OB 

(4) lim (u(x)) = 0} 

x— >oo 

equipped with the bundle projection r(b,b,u) = b and the vector bundle structure 

(b, b\,ui) + (b, b 2 , u 2 ) = (b, b\ + b 2 ,u 1 + u 2 ) 

for all (b,bi,Ui),(b,b 2 ,u 2 ) G r~ l {b) and all b G Emb(B). Additionally, the map 
p : (TQ/G) ->■ TEmb(£>) given by p(b,b,u) = (b,b) is well defined. 

Proof. Let us temporarily label the set in question by U A" and define the map 
tt /g : TQ -> A by 

tt/g(6, b, f, ip) = (b,b\ipof- 1 ). 

We see that tt/g(v o t/;) = ir/ G (v) for all G G and v G TQ. Therefore, ir / G 
maps the coset v ■ G to a single element of A. Conversely, given an element 
(b, b,u) G A we see that tt^(6, b, u) is the set of element in (6, b, ip, ip) G TQ such 
that u = ipoip^ 1 . However, this set of elements is just the coset v-G where v is any 
element such that tt/g(v) = (b,b,u). Thus tt/g induces an isomorphism between 
TQ/G and A. Additionally we can check that ir/ G (v + w) = n/ G (v) + n/ G (w) and 
t { 71 /g{ v )) — t q{ v ) ' G. Therefore, the desired vector bundle structure is inherited 
by TQ/G as well and tx/ g becomes a vector bundle morphism. Finally, the map 
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p(b, b,u) = (b, b) is merely the map T(3 : TQ — > T(Q/G) divided by G. That is to 
say poj[ G = T/3. This equation makes p well defined because (3 is C7-invariant. □ 

Just to reiterate. The fluid velocity component u in (b, b, u) is not a vector-field 
(strictly speaking) because the boundaries of u may point in directions outside of 
it's domain, r^. This reflects the fact that the domain of u is time dependent. 

We now define the reduced Lagrangian I : TQ/G — > M by 

l(b,b,u) = L B (b,b) + ||n(a;)|| 2 rfvol(a;). 

This induces the standard Lagrangian L \= I o-k j G : TQ — > R which is G-invariant 
by construction. Additionally we wish to add a viscous force on the fluid, : 
TQ T*Q. We first define the reduced force field /„ : TQ/G (TQ/G)* by 



(fv(b,v b ,u),(b,w b ,w)) = fi / (Au(x),w(x))dvol(x), 

and define : TQ(b ) -+ T*Q(b ) by 

(F^(v),w) = (f^n /G (v)),7T /G (w)). 

We finally define the total force on our system to be F = F^ + (Fg o T/3) where Fq 
is the dissipative force on the shape of the body mentioned in §3.2| This total force 
F descends via vr /G to a reduced force F /G : TQ/G ->■ (TQ/G)* where (TQ/G)* 
is the dual vector bundle to TQ/G. The reduced force is given explicitly in terms 
of f p and F B by F/ G = / M + (F B o p). One can verify directly from this expression 
that (F(v),w) = (F /G (7r /G (v)), 7i /G (w)). 

We now introduce a consequence which follows from G symmetry of F and L. 

Proposition 4. Let Flow* : TQ — > TQ denote the time-t evolution operator 
corresponding to the Lagrange- d'Alembert equations induced by the Lagrangian 
L : TQ — > R and the force F : TQ — > T*Q. Then there exists a unique evo- 
lution operator flow* : TQ/G — » TQ/G such that flow* o tt/ g = tt/ g o Flow*. In 
other words, flow* is the unique evolution operator which makes the diagram 

T Q ^ TQ 



TQ/G TQ/G 



commute. 



Proof. Let q : [0, t] — > Q be a curve such that the time derivative (q, q) : [0, t] — > TQ 
is an integral curve of the Lagrange-d'Alembert equations with initial condition 
(q,q)(0) and final condition (q,q)(t). Then the Lagrange-d'Alembert variational 
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principle states that 

5 f L((q,q)(r))dT= f (F((q,q)(rl5q(r))dr 
Jo Jo 

for all variations of the curve g(-) with fixed endpoints. Note that for each if) G G 

the action satisfies f* L((q,q)(r)dT = f Q L((q, q)(r) oip)dr and the variation on 

the right hand side of the Lagrange-d'Alembert principle is 

1 (F((q, q)(r)), 5q(r))dr = [ (F /G ( VG ((q, q)(r), 7r /G (5q (r)))dr 
o Jo 

= [ (F((q,q)(r)oif}),8q(r)oif;)dT. 
Jo 

Therefore we observe that 

5 f L((q,q)(r)o^)dT= [ (F((q, q)(r) o if}), 5q(r) o ^)dr 
Jo Jo 

for arbitrary variations of the curve q(-) with fixed end points. However, the 
variation Sq o ip is merely a variation of the curve q o ^(-) becuase 

d 



5q(r) o ip 



(g(r, e) o ijj) 



e=0 

and if q(r, e) is a deformation of q(r) then q(r, e) o ijj is a deformation of g(r) o ip 
by construction. Therefore, 

'* L((q, q)(r) o ^)dr = f (F((q, q)(r) ° ^), 5{q o ^)}dr 

Jo 

for arbitrary variations of the curve q o ip with fixed end points. This last equation 
states that the curve (q, q) o ip satisfies the Lagrange-d'Alembert principle. Thus 
Flow*((g, q)(0) o if}) = (q,q)(t) o ip. Since if) G G was chosen arbitrarily we may 
apply Flow* to the entire coset (q, q) • G to find 

(5) Flow'((g, q) ■ G) = Flow*(g, q) ■ G. 

This last equation describes a map from cosets to cosets, that is from TQ/G to 
TQ/G. In other words, equation ^ is the map flow* : TQ/G — >■ TQ/G. □ 

3.4. Reduction by frame invariance. Consider the group SE(d) consisting of 
rotations and translations of lR d . Each z G SE(<f) sends (b,b,u) G TQ/G to 
(z(b,b), z*u) G TQ/G where is the push-forward of the fluid velocity field to 
the domain pza(z-b). This action is free and proper on TQ/G so that the projection 
[•] : TQ/G -> [TQ/G] where 

'^1 - i| 
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is a principal bundle [A MQ O, Prop 4.1.23]. Additionally, SE(d) acts by vector 
bundle morphisms which are isomorphisms on each fiber so that [TQ/G] inherits 
the vector-bundle structure from TQ/G. 

Proposition 5. There exists a function [I] : [TQ/G] — > R such that I o [■] = [I], 
a vector-bundle projection [r] : [TQ/G] — > [Emb(B)], a map [p] : [TQ/G] — > 
T[Emb(i3)] such that the diagrams 

TQ/G [TQ/G] TQ/G -^U [TQ/G] 

[p] 

Emb(£) — [Emb(B)] TEmb(B) -^L> T[Emb(£)] 

commute. Lastly, there exists a force [F/q\ '■ [TQ/G] — > [TQ/G]* define by the 
equation ([F/g\([€\), [v]) = (F/g{Q,v) f or an U £iV e TQ/G above the same base 
point. 

Proof. We first show the existence of [I]. Let z G SE(d) and £ = (b,b,u) G TQ/G. 
Note that when z acts on a vector (x, 5x) G TIR d it will preserve the length of the 
the vector. That is to say ||5x|| 2 = ||z ■ 5x\\ 2 . Thus we see that 

l{z ■ = L B (z ■ {b, &)) + [ \\z,u{x)\\ 2 d vol 
= L B (b,b)+ I |Kx)|| 2 dvol 

= K0 

where we have used the change of variables formula and the fact that L B is SE(d) 
invariant by assumption. Since z is aribtrary we find that applying I to the entire 
coset SE(cf) • £ yields a single real number. This means we may define a function 
[/] : [TQ/G] — > K by the condition I — [I] o [■]. Continuing in this manner we see 
that 7"(z-£) = z-t(£) so that there is a well defined map [r] : [TQ/G] — > [Emb(B)] 
defined by the condition r = [r] o [•] and the analogous argument map be applied 
to the map p : TQ/G -> TEmb(£) to yield a map [p] : [TQ/G] ->■ T[Emb(<B)]. 
Lastly, we note that 

(F /G (* • £), (2; • 77)) = (F B (z -(b,b),z- (6, 6') + / (A(z*u)(x), z*u'(x))dvol 

JXXzb 



= (F B (b,b),(b,b')) + J (Au(x),u'(x))dvol 
= (F /G (0,V)- 

Where £ = (b, b, u) and r\ = (b, b', u') are arbitrary element of TQ/G over the same 
fiber and z is an arbitrary element of SE(d). Since z is arbitrary, there must exist a 
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unique function [F /G ] : [TQ/G] [TQ/G]* such that ([F /G ]([£]), [rj]) = (F /G (0,v) 
for arbitrary £, G TQ/G. □ 

In the preceding proof we also illustrated a few other things. In particular, the 
Lagrangian L and the unreduced force F are also SE(cf) invariant under the left 
action z • (b, ip) = (z • b, z*ip) for each z G SE(d) and (6, ip) G Q. Therefore the 
flow of the system is also invariant with respect to the action of SE(eZ) and so we 
may define a flow on a space which has been quotiented a second time by this left 
action. 

Proposition 6. Let flow* : TQ/G — > TQ/G be the time t evolution operator of 
Proposition^ There exists an evolution operator [flow]* : [TQ/G] — > [TQ/G] such 
that the following diagram 

TQ/G ^> TQ/G 



[TQ/G] [ T Q/ G] 



commutes. 



Proof. Let q : [0, T] — > Q be a curve such that the time derivative (q, q) : [0,T] — > 
TQ is an integral curve of the Lagrange-d'Alembert equations with initial condi- 
tion (q, q)o and final condition (q, q)x- Then the Lagrange-d'Alembert variational 
principle states that 



6 [ L(q,q)dt= [ (F(q,q),5q)dt 
Jo Jo 







for all variations of the curve q(-) with fixed endpoints. Note that for each ip G G 
and z G SE(d) the action satisfies L(q,q)dt = J Q T Liz ■ (q, q) o ip)dt and the 
variation in the work is 

(F(q,q),6q)dt= \ (F /G (z ■ TT/ G (q, q), z ■ ir/ G (6q))dt 

T 

(F(z ■ (q, q) o if)) 7 z ■ 8q o %fj)dt. 

o 

Therefore we observe that 

f T f T 

5 / L(z ■ (q,q) o ipjdt — / \F(z ■ (q,q) ° ip), z ■ Sq o ijj)dt 

Jo Jo 
for arbitrary variations of the curve q(-) with fixed end points. However, the 
variation z ■ Sq o ip is merely a variation of the curve z ■ q o ip(-). Therefore, 

6 f L(z- (q, q) o ip)dt = [ (F(z ■ (q, q) o ip), 5(z ■ q o ip))dt 
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for arbitrary variations of the curve z-qoip with fixed end points. This last equation 
states that the curve z-(q,q)oif) satisfies the Lagrange-d'Alembert principle. Thus 
Flow T (2; • (q, q)o o if)) = z ■ (q, q)x ° if)- Since if) G G was arbitrary we may apply 
Flow T to the entire coset z ■ (q, q) - G to find 

F\ow T (z ■(q,q)-G) = z- Flow T (g, q) ■ G. 

This last equation states that 

flow T (z ■ i) = z ■ flow T (0 

where £ = vr/ G ((g, q) ). However, as z G SE(d) is arbitrary we can multiply by all 
of SE(d) to get 

(6) flow'(SE(rf) ■ = SE(d) • flaw'(f). 

This last equation describes a map from cosets to cosets, that is from [TQ/G] to 
[TQ/G]. In other words, equation @ is the map [flow]* : [TQ/G] [TQ/G]. □ 

4. Asymptotic Behavior 

Our experience of observing a dumpling floating in a bowl of soup should sug- 
gests that the passive motion of the system tends towards one where the fluid 
is stagnant and the shape of the solids has settled to some minima of the elastic 
potential energy. In this section we will confirm this using language presented thus 
far. 

Proposition 7. Let q : [0, oo) — )■ Q be a curve such that the time derivative 
(q,q) : [0, oo) — > TQ is an integral curve of the Lagrange-d'Alembert equations for 
the Lagrangian L and the force F. Then the u -limit set of (q, q)(-) is contained in 
the set 

dU- 1 (0):={(q,0)eTQ | dU(q) = 0}. 

Proof. Roughly speaking the proof goes as follows. We find that the energy is a 
quadratic positive definite funtion. The time derivative of the energy is negative as 
along as the velocity is non-zero. This is due to dissapation from viscous frictions. 
The positive-definitness of the energy combined with its negative time derivative 
will allow us to conclude that the velocity is sent to zero. That is to say, the system 
evolves towards the zero section of TQ. We will then find that acceleration on 
the zero- section vanishes if and only if dU. This means that the system goes to 
points in the zero section where dU vanishes and concludes the proof. 

We will now proceed to prove this in a more explicit manner. Before we begin 
we will introduce some handy notation. Given a vector bundle tt : V — > M there 
exists a unique section called the zero-section which maps each m G M to the 0- 
vector in V above the point m. We will denote this section for an arbitrary vector 
bundle by (-)2. This point is minor perhaps, becase it is customary to identify the 
zero-section of a vector bundle with the base of the bundle. However, to avoid 
causing confusion we will not make this identification. 
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The energy is the function E : TQ — > R given by 
E(q,q) :={WL(q,q),q)-L(q,q) 

= \( p B {x)\\b{x)fdvo\ B + \ f ||0(x)|| 2 dvol - U(b) 

Given any Lagrangian system on a Riemannian manifold where the Lagrangian is 
the kinetic energy minus the potential energy, the time derivative of the generalized 
energy under the evolution of the Lagrange-d'Alembert equations is given by E = 
(F(q), q). In this case we find 

E(q,q) = (F B (b),b) + (F,(q,q),(q,q)). 

However, the right hand side is a (weakly) positive definite quadratic form on each 
fiber of TQ. Therefore the w-limit of (q, q)(-), denoted must be a subset of the 
zero section of TQ, which is identifiable with Q itself. In other-words, M w C Q. 
Let (6, ip) = q G M w (g, q). Then the Lagrange-D'Alembert equations is equivalent 
to 

^ = W(&) + tt(FM)) 

where Jj : T*Q — > TQ is the sharp map associated the metric on Q. However, 
F{q,q) = when q = 0, which is the case for point in M w (g,g). Thus, the 
evolution of the system is given by the eqaution 

for points in the zero section of TQ. The covariant derivative takes place above a 
vector with zero velocity so that q G TTQ is a vertical vector. The above equation 
says that the vertical part of q is given by VU. There the vector field on TQ is 
transverse to the zero section of when dU ^ 0. Hence we find that we can restrict 
M u further. That is to say, M u C dU~ l {Q). □ 

Corrollary 1. Let [U] : [Emb(B)] — > 1R be the unique function on the shape- 
space of the body such that [U]([b]) = U(b) for all b G Emb(i3). Assume that 
[U] has a unique minimizer s m j n G [Emb(i3)]. Then if (q,q) : [0, oo) — > TQ is 
an integral curve of the Lagrange-d'Alembert equations, then [£](t) = [vr/c^, <l{t))} 
must approach (s m i n )° G [TQ/G]. If the flow of the Lagrange-d'Alembert equations 
is complete, this means that (s m i n )° is a global (weakly) hyperbolically stable fixed 
point for the vector field X\tq/g\ ■ 

Proof. In proposition [7] we showed that solutions approach points within the set 
dU~ l (Qi) asymptotically. This implies that the dynamics on [TQ/G] must approach 
g?[£7] -1 (0) asymptotically. However, there is only one such point. □ 

We can interpret corrollary [T] to mean that "a motionless body in stagnant water 
is a stable state for dead fish". In the next section we will periodically perturb 
this stable equilibria to obtain a loop in [TQ/G]. 
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5. Swimming 

When one activates one's muscles they do so by changing chemical potentials 
which then stiffen muscle fibers, causing them to contract. In terms of Lagrangian 
mechanics, this has the effect of altering the elastic potential energy of the muscle 
tissue. Therefore to model the periodic muscle activation of a fish swimming we 
may consider a time periodic potential energy U : S 1 x [Emb(S)] — > 1R, which we 
define on shape space. If we add this potential energy to our Lagrangian we do 
not break any of the symmetries discussed thus far. However, this does alter the 
dynamics. In particular, adding U to the Lagrangian is equivalent to adding a 
time periodic force F = dU to the equations of motion. If we let X T q 6 %(TQ) 
denote the vector field for the flow, Flow', then the flow in the time-periodic 
augmented phase space is given by the vector field (d t ,X T g) 6 X(5' 1 x TQ). 
Adding a periodic potential energy effects the dynamics in augmented phase space 
by the addition of a vector-field Y TQ : S 1 x TQ ->■ TTQ |MR99l §7.8] . In particular, 
the dynamics of the periodically perturbed system in augmented phase space are 
given by (d t ,X T Q + Y T q) e X(S' 1 x TQ). If U only depends on the shape of 
the body, then Ytq is both G and SE(ci) invariant and thre exists a vector field 
Y[tq/g\ ^ £([TQ/G]) so that the dynamics on [TQ/G] in augmented phase space 
are given by (d t ,Y [TQ/G] + X [TQ/G] ) e X(5' 1 x [TQ/G]). Because X [TQ/G] exhibits 
an asymptotically stable fixed point , x G [TQ/G], then the orbit {(x ,t) : t e 
S 1 } C S 1 x [TQ/G] must be a (weak) hyperbolic stable limit cycle for the vector 
field (d t , X[ T q/g]). If we were observing a finite dimensional manifold we could 
invoke the persistence theorem to state the existence of a perturbed limit cycle in 
the dynamics of (d t , X\tq/g] +Y\tq/g\) for sufficiently small oscillations. However, 
this final leap is problematic. 

5.1. Analytical Issues and a work-around. Unfortunately, our manifold is 
infinite dimensional and Frechet. I am unable to assert anything with absolute 
certainty on the existence of limit cycles in [TQ/G] for arbitrary perturbations. 
Specifically, the strength of the stable fixed point of xq is tied to the spectrum 
of the Laplacian operator on r^, which is strictly negative but approaches 0. In 
the language of |Fen72j . there is no spectral gap condition. Secondly, even if there 
were a spectral gap, the existence of a limit cycle involves fixed point theorems 
which assume that our space is complete. 

However, perhaps these problems are not so devastating. There exists a number 
of finite dimensional models for the space TQ/G used by engineers to study fluid 
structure interaction. It is fairly common to approximate the fluid velocity field on 
a finite dimensional space and model the solid using a finite element method. For 
example consider the immersed boundary method |Pes02] . Let us call this space 
(TQ/G) discrete- Moreover, usually one can act on (TQ/C7)discrete by SE(<i) simply 
by rotating and translating the tetrahedra in the case of finite element methods, or 
rotating the basis functions for a spectral method. If the model on {TQ/G) discrete 
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converges as the time step goes to zero then we could reasonably restrict ourselves 
to methods which dissipate energy at a rate which is quadratic and positive definite 
in velocity. This is not too much to expect because a good method ought to 
converge and [Pes02j or [GM77J are both methods which exhibit this property. By 
the same arguments as before the dynamics will exhibit a hyperbolically stable 
equilibria on the quotient space ^ T ^^|^ scrctc . Upon adding a periodic perturbation 

to the dynamics on ^ TQ ^|^ scrctc one could apply the persistence theorem directly 
to assert the existence of a hyperbolically stable limit cycle in augmented phase 
space. If the model on ^ T< ^^|^ scr ° tc converges as the grid size goes to zero, then for 
a finite grid resolution there is a trajectory in TQ/G which is well approximated 
by the definition of the term "convergence". Therefore one could say that this 
limit cycle approximates the real dynamics, although perhaps only on a finite time 
span. That is to say, there is a loop in [TQ/G] which nearly satisfies the exact 
dynamics. Therefore instead of searching for conditions under which a limit cycle 
exists for the exact dynamics we shall explore the consequence of evolving on loops 
in [TQ/G] which approximate the exact dynamics. Given one of these loops we 
will be able to interpret regular motion (a.k.a swimming) by the use of phase 
reconstruction formulas induced by lifting the loop in [TQ/G] to a path in TQ/G. 



t=2 




B 



t=o t=i 



t=2 



t=3 



Figure 4. The bottom row (B) depicts snapshots of a periodic 
trajectory in [TQ/G]. As the trajectory is periodic, all the snapshots 
are identical. The top row (A) is a lift of the trajectory in (B) 
to TQ/G. The top row appears to be steady swimmming. Each 
snapshot is related to the previous by a (fixed) rigid rotation and 
translation. 

5.2. Swimming as a limit cycle. Let [7] : S 1 — > [TQ/G] be a loop such that 
graph([7]) := {{t, [y](t)) : t £ S 1 } is approximates the dynamics in the augmented 
phase space S 1 x [TQ/G]. There must exists a curve 7^ : [0,1] — > TQ/G such 
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that [7 (1) (t)] = [7KO for t G [0,1]. Additionally, since [7](0) = py](l) there must 
exists a z G SE(<2) such that 7^(1) = z ■ 7^(0). However, we see that z ■ 7^(t) 
is a smooth curve which projects to the loop [7]. Therefore, by concatenating the 
curves 7^(-) and z ■ 7^(-) we can get a longer integral curve ^ 2 > : [0, 2] — > TQ/G 
defined by 

^2 )(t) = {i (1) (t) for te [0,1] 

7 U \z-^(t-l) forte (1,2]' 

Note that 7 (2) (2) = z ■ 7 (1) (1) = z 2 ■ 7 (1) (0). Observing that Z7 (2) (-) projects to [7] 
as well we can extend 7^(-) to get a curve 7® : [0, 3] — > TQ/G. By induction we 
can extend this argument indefinetly to get a curve for all positive time 

7 (t)=* LtJ -7 (1) (*-Lt|) 

where [ij = sup{fc eZ:Kt}. One way to interpret this visually is as follows. If 
one were to take a snap shot of the body immersed in the fluid every second, then 
one would observe the same picture over and over again, except each picture would 
be translated and rotate by the fixed element z G SE(d) (see figure [4p. Moreover, 
since [7](-) is a stable limit cycle, this behavior is stable and so we should expect 
to observe it in real systems (see figure [I]). Of course, it is possible that the 
reconstructed phase-shift z G SE(d) is simply the identity, which would mean 
there has been no net motion after one flap of the fins. The main point is that this 
regular behavior is stable and generally z will not be the identity. Constructing 
ways of controlling the phase-shift z will certainly bring us into future work, which 
is the topic of the next section. 

6. Conclusion and Future work 

It is widely observed that steady swimming is periodic, and this observation 
inspired the question "is it possible to interpret swimming as a limit cycle?" . In this 
paper we have answered this question in the affirmative. We did this by carefully 
constructing the phase space and then reducing it by the appropriate symmetries. 
We arrived at a system where we could invoke the persistence theorem to construct 
loops that approximate the dynamics in the reduced phase space. Lifting these 
loops by phase reconstruction formulas allowed us to assert the existence of stable 
trajectories which strongly resemble swimming. The stability and the regularity of 
these trajectories were important features, as both of these features are observed 
in real systems. 

Given the complexity of fluid-structure interaction it is not immediately clear 
that one should have expected such orderly behavior. This orderliness has the 
potential to be exploited in a number of applications. 

(1) Robotics and Optimal Control The interpreation of swimming as a 
limit cycles may permit a new paradigm for controller design. For example, 
if we assume that we have a potential energy which is a sum of two parts 
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£/([&]) + f/ c ([6]) where the U c term depends on a control parameter c G C. 
We can denote the set of loops in C by loop(C) then we may consider the 
map T : loop(C) — > loop([TQ/G]) which outputs the periodic limit cycle 
in [TQ/G] which results from using a loop in C. We may use T to define 
our cost function on the space loop(C). Such a cost functional would be 
one which is sensitive to the long term behavior and does not over-react to 
the transient dynamics. 

(2) Transient dynamics Even though trajectories may approach a limit cy- 
cle, the transient dynamics will matter. The transient dynamics will re- 
orient and translate the body before any orderly periodic behavior kicks 
in. Therefore, if one desires to create locomotion through periodic control 
inputs, one should try to get onto a limit cycle quickly in order to min- 
imize the duration where the complicated transient dynamics dominate. 
Both these areas can be studied if and only if the map T exists. A main 
contribution of this paper a demonstration of an "approximate" map T. 

(3) Pumping In the current setup one could consider a reference frame at- 
tached to the body. In this reference frame "swimming" manifests as fluid 
moving around the body in a regular fashion. This change in our frame 
in reference is really a description of pumping and suggests the stability of 
fluid pump may be proven in a similar manner. 

(4) Passive Dynamics This paper does not address the dual problem. By 
the dual problem we mean: "Given a constant fluid velocity at infinity, 
what periodic motion (if any) will a tethered body approach as time goes 
to infinity?" In this dual problem, the motion of the body is given first, 
and parameters such as the period of the limit cycle are emergent phe- 
nomena. In particular, the dual problem of a flapping flag immersed in a 
fluid with a constant velocity at infinity has received much attention in the 
applied mathematics community (see [SVZ05J and references therein). It 
is aboserved that stable periodic motion occurs for a range of velocities at 
infinity while resonant modes of the flag lead to chaotic behavior in certain 
regions. This problem has been studied from a variety of angles already. 
However, it may be insightful to view it in a manner similar to the analysis 
presented in this paper. Such an interpretation could help generalize well 
known results which would otherwise be specific to flapping flags. 

(5) Other types of locomotion The notion that walking may be viewed 
as a limit cycle has been around for a while (see [HW07J and references 
therein). Moreover, it is conceivable that flapping flight is a limit cycle 
as well. However, for both these systems SE(3) symmetry is broken by 
the direction of gravity. Because of this, it is not immediately clear that 
one can import the methods used here to understand flapping flight and 
terrestrial locomotion. However, perhaps this is merely a challenge to be 
overcome. In particular, these systems still exhibit SE(2) symmetry. So it 
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is conceivable that there exists a fixed point for a reduced system which 
is hyperbolic but with an unstable mode. Using the same construction 
one can infer the existence of a periodic orbit which is not a limit cycle. 
It would be interesting to see if we could describe the unstable directions 
physically and inspire controls based upon these insights. 
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